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Abstract 

We study the p-meson spectral function in hot nuclear matter by taking into account the isospin- 
symmetric pion and the nucleoli loops within the quantum hadrodynamics (QHD) model as well 
as using an effective chiral SU(3) model. The spectral function of the p meson is studied in the 
mean field approximation (MFA) as well as in the relativistic Hartree (RHA) approximation. The 
inclusion of the nucleon loop considerably changes the p-meson spectral function. Due to a larger 
mass drop of p meson in the RHA, it is seen that the spectral function shifts towards the low 
invariant mass region, whereas in the MFA the spectral function is seen to be slightly shifted 
towards the high mass region. Moreover, while the spectral function is observed to be sharper with 
the nucleon-antinucleon polarization in RHA, the spectral function is seen to be broader in the 
MFA. 
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I. INTRODUCTION 



The study of strongly interacting nuclear matter under extreme conditions has attracted 
a lot of attention during the recent years, both theoretically and experimentally. The prop- 
erties of hadrons at high temperatures and densities are quite different from the properties 
of hadrons in vacuum. The quantum hadrodynamics QHD-I model (Walecka model) and 
its extensions have been widely used to discuss the properties of the nuclear matter and 
finite nuclei [1-5]. The ongoing relativistic heavy ion collision experiments, at the high en- 
ergy accelerators SPS, CERN, Switzerland; SIS, GSI, Germany; RHIC, BNL, USA; LHC, 
CERN, Switzerland, and the compressed baryonic matter (CBM) experiments planned at 
the future facilities at GSI, Germany, are intended to probe matter at high temperatures 
and densities. The hadrons modified in the hot and dense hadronic medium resulting from 
heavy ion collision experiments, affect the experimental observables. E.g., the dilepton spec- 
tra observed from heavy ion collision experiments at the SPS [6, 7] are attributed to the 
medium modifications of vector mesons [8-15], and can not be explained by vacuum hadronic 
properties. 

It is predicted that in ultra-relativistic heavy-ion collisions quark-gluon plasma may be 
produced and the spontaneously broken chiral symmetry may be restored. Amongst the 
proposed signals for detecting the quark-hadron phase transition, dileptons and photons are 
considered to be the cleanest ones because they do not interact with the hadronic medium 
and emerge from the heavy ion collision experiments almost undisturbed [16-18]. Also 
considering the fact that the light vector mesons can directly decay to dilepton pairs, the 
study of the p meson in the medium is interesting because of its relatively large decay width 
as compared to those of u and mesons. [19-25]. 

One of the hot topics in recent years is discussing the low invariant mass dilepton produc- 
tion in heavy-ion collisions. The dilepton distribution is related with the spectral function of 
p meson [22, 26]. In the present work, the p meson spectral function is calculated accounting 
for effects of the pion as well as the nucleon loops. The effect of modification of the nucleon 
mass on the p meson spectral function is discussed within the framework of both the Walecka 
model and a chiral SU(3) model. In both the models, the effective nucleon mass is modified 
with density and temperature and this in turn modifies the p meson spectral function. 
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We organize the paper as follows: Section II gives a brief description of the QHD-I 
(Walecka) model and the chiral SU(3) model used in the present investigation. In Section 
III, we discuss the nucleon properties in hot and dense matter in both of these models. 
Section IV discusses the effects of finite temperature and density on the spectral function of 
the p meson. Section V contains the results and discussion and in section VI, we summarize 
our main findings of the present investigation and discuss possible outlook. 

II. THE HADRONIC MODELS 
A. QHD-I (Walecka) model 

The Lagrangian density for the model is given by [4, 5] 

L= i> [ 7 „ (id* - — (M — g a a)\ ^ + \ (d^&a - m 2 a a 2 ) - |/*7 3 - ^Xa 4 

(1) 

where ip is the nucleon field, a is the neutral scalar meson field and oj is the isoscalar vec- 
tor field. F^ v = d^oj" — dyoo^ is the field tensor for the vector meson, oj and SL contains 
counterterms used for renormalization. The parameters M, g a ,g u , m a , m u , k, and A are phe- 
nomenological constants that are determined from the nuclear matter saturation properties. 

In the QHD-I (Walecka model), the nucleons interact through the exchange of a and oj 
mesons. The a exchange gives the attractive force while the u exchange attributes to the 
repulsive interaction between the nucleons. We use the QHD-I model to obtain the effective 
nucleon mass and effective chemical potential p* in the hot and dense nuclear matter 
and use the obtained and p* in mean field or relativistic Hartree approximation for 
investigating the in-medium properties of the vector mesons. 
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B. The hadronic chiral SU(3) x SU(3) model 



The effective hadronic chiral Lagrangian density used in the present work is given as 

L = Lkin + Lew + L vec + L + Lsb (2) 

W=X,Y,V,A,u 

Equation (2) corresponds to a relativistic model of baryons and mesons adopting a nonlinear 
realization of chiral symmetry [27-29] and broken scale invariance as a description of the 
hadronic matter. Here, Lj, in is kinetic energy term, Lbw is the baryon- meson interaction 
term in which the baryons-spin-0 meson interaction term generates the baryon masses. L vec 
describes the dynamical mass generation of the vector mesons via couplings to the scalar 
mesons and contains additionally quartic self-interactions of the vector fields. L contains 
the meson-meson interaction terms as well as a scale invariance breaking logarthimic po- 
tential. L S b describes the explicit chiral symmetry breaking. The baryon-scalar meson 
interactions generate the baryon masses and the parameters corresponding to these inter- 
actions are adjusted so as to obtain the baryon masses as their experimentally measured 
vacuum values. For the baryon-vector meson interaction terms, there exist the F-type (an- 
tisymmetric) and D-type (symmetric) couplings. Here we use the antisymmetric coupling 
[30-32] because, following the universality principle [34] and the vector meson dominance 
model, one can conclude that the symmetric coupling should be small. Additionally we 
choose the parameters [30, 33] so as to decouple the strange vector field M ~ sj^s from the 
nucleon, corresponding to an ideal mixing between u and <fi. A small deviation of the mixing 
angle from the ideal mixing [35-37] has not been taken into account in the present inves- 
tigation. The Lagrangian densities corresponding to the interaction for the vector meson, 
L vec , the meson-meson interaction L and that corresponding to the explicit chiral symmetry 
breaking L SB have been described in detail in references [30, 33]. 

To investigate the hadronic properties in the medium, we write the Lagrangian density 
within the chiral SU(3) model and determine the expectation values of the meson fields by 
solving the equations of motion of the scalar fields at finite temperature and density obtained 
by minimizing the thermodynamic potential. 
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III. NUCLEON PROPERTIES IN HOT NUCLEAR MATTER 



First we proceed to study the hadronic properties in the QHD-I model. We study the 
nucleon properties in both the mean field and relativistic Hartree approximations (RHA). 
In the mean field approximation (MFA), the meson field operators can be approximated by 
their expectation values, which are the classical fields. The expectation value of the scalar 
field, <T shifts the nucleon mass from M N to = M N — g a a and the vector field gives 

2 2 

rise to an effective chemical potential, pf — p — 9w £ pB ■ The nucleon effective mass can be 
determined self-consistently by [4] 

M N -M* N = -^p s B (3) 
where the scalar density p s B of the nuclear matter is 

a£ = 7"/^||(Mp)+Mp)) (4) 

with the spin-isospin degeneracy factor, 7 = 4 for symmetric nuclear matter. The effective 
chemical potential p* is determined by 

2 

p = p — (5) 



m 2 v 



with p B as the baryon density, given by, 

/d^p 
j^(n N (p) - n N (p)) (6) 

Here, and n N are the thermal distribution functions given as 

UN ^ = e (E* N -X T + V ™M = e (E* N+ X T + 1 - ( 7 ) 



where, E* N = + (M* N ) 2 is the single particle energy of the nucleon and T is the tempera- 
ture. One can solve the above coupled equations numerically to obtain the effective nucleon 
mass and the effective chemical potential p* for given density, pb and temperature, T. 
The RHA takes into account the vacuum fluctuation corrections to the mean field results. 

Next we proceed to study the hadronic properties in the chiral SU(3) model used in the 
present investigation [38, 39]. The Lagrangian density in the mean field approximation is 
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given as 



L BX + L BV = -ip N [g Nw ~iQU + M* N ] ip N , 

t X , ,2 i „4, ,4 



^ Xo 

L = \hx 2 (« 2 + C 2 ) - h (a 2 + c 2 ) 2 - h + C A 



hx^ 2 C 



+k,X A + \x A ^ X 



Xt 3 



X In 



Lsb — \ — 

\Xo, 



m 



lf n a + I V2m 2 K f K - -^ m lfn) C 



(8) 



where = —QaN ' — 9cnC is the effective mass of the nucleoli. The thermodynamical 
potential of the grand canonical ensemble Q per unit volume V at a given chemical potential 
ji and temperature T can be written as 

n r d 3 P 



y = —L vec — L — L S b — v vac + In J 
d 3 p 



(27T) 



■ E N(p)( n N(p) +n N (p)) 



-7iv 



(27T) 



-u* N (n N (p) - n N {p)) 



(9) 



Here 7^ are the spin-isospin degeneracy factor, and 7^ = 4 for symmetric nuclear matter. 
The tin and njy are the thermal distribution functions for the nucleon and the antinucleon 
given in terms of the effective single particle energy, E* N) and the effective chemical potential, 
/j,*, as given by equation (7). The mesonic field equations are determined by minimizing the 
thermodynamic potential. We shall use the frozen glueball approximation (x = Xo) > since 
the dilaton field which simulates the gluon condensate changes very little in the medium. 
We then have coupled equations only for the fields a, ( and u as given by 



9 (n/V) = foxV - Ah (a 2 + C 2 ) o - 2k 2 a 3 - 2k sX a( - 2*** 



da 



d(n/v) 



3a 



9 dM* q 
+m 2 J n + -^-p% = 0, 

= A; oX 2 C-4A; 1 (a 2 + C 2 )C-4A; 2 C 3 -A;3xa 2 

1 



(10) 



$X*_ 
3C 



V2m 2 K f K - -^ m lU 



= 0, 



djsyv) 

du; 



= -m z w uj - Ag\uj 3 + g Nu p B = 0, 



(11) 



(12) 



which have to be solved self-consistently to obtain the values of a, ( and u . Here p s B and 
Pb are the scalar and vector densities for the nuclear matter at finite temperature, T given 
by equations (4) and (6). 



IV. RHO MESON SPECTRAL FUNCTION 



The most consistent approach for studying the hadronic matter produced in ultra- 
relativistic heavy-ion collisions at high temperature and density is the finite temperature 
field theory [40, 41]. In the present calculation thermal effects enter through thermal nu- 
cleon and pion loops. In Minkowski space, the self-energy of the p vector meson can be 
expressed as [22, 42] 

IP" (A;) = n L (k) If + U T (k) PF, (13) 

where k 2 = k 2 ,— \k 2 \. The P£ v and P^ v are the longitudinal and transverse projection tensors 
defined as 

pQ° = p°* = = 0, P« = 5 ij -ih, P? = - - P% v . (14) 

\k 2 \ 

Ii L and Il r are related to the components of the self-energy by 

U L (k) = ~^(k), 

Ht (k) = l(K~ U l (*)) • (15) 

The imaginary part of the retarded propagator is referred to as the spectral function 
and is related to the dilepton production. The study of the p meson spectral function is 
attributed to calculating the in-medium self-energy of the p meson. 



A. pNN interaction 

The contribution of nucleon excitations through nucleon-loop to p self-energy is analyzed 
in terms of the effective Lagrangian density [43, 44] 

(16) 
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where V£ is the p meson field and \I/ is the nucleon field. The one-loop contribution to the 
polarization tensor is 



r M (*) 



1 



7^ - M* N 



r„ (-k) 



1 



r (p-k) u -M* N 



where 



r» (k) = ^ + 



(17) 



(18) 



with a^ u = | [7^, 7^], M^r and are the nucleon masses in vacuum and in the hot hadronic 
medium respectively. The polarization tensor (k) can be separated into two parts, 

K N (*) = (^r " ^) n ^ (*) + n C (*) > (19) 

corresponding to the vacuum and the matter contributions. Using dimensional regularization 
and taking a phenomenological subtraction procedure [44], the vacuum part (T = 0) is, 



where 



Jo Jo 



dx In C, C = „ 

M 2 , - x (1 - x) A; 2 



(21) 



The second part n^^(A;) is given in terms of the thermal distribution functions ([/,*, T) 
and un (fi*,T). The longitudinal and the transverse parts of the self-energy can be calculated 
from Ii^ v using equation (15). The longitudinal part of the matter part of the self-energy is 
given as: 



Kl (*) = n?£2 (*) + nSS2 (*) + n&J (*) , 



T pATiV 



where, 



(22) 



'k 2 + 4(E%) 2 ' 



p\k\ 



n£S2(fc) 



2M, 



f 

Jo 



pdp 



E 



n N (//, T) + njy (/A T)) In A, 



A? 



( 9 2 p nn \ ( k 2 
I 4tt 2 i 



8fc 2 - 



V|A;| 2 

h 2 h 2 



p 2 dp 
~E~ 



n N (n*,T) + uat (/i*,T)) 



A' 



In A 



4K* A- 2 A- 2 



(23) 



Here p — \p\ and E* N = \jp 2 + Mff. k is the 3- momentum of the vector meson, p. A and 
B are defined as 

. ->.\ 2 



A 



B = 



(k 2 + 2p\k\) -A{E* N ) 2 k 2 
(k 2 - 2p\k\) 2 - A(E* N yk 2 
k*-A(p\k\ + E* N k { 
k*-4(p\k\-E* N k ) 



(24) 



Our longitudinal results agree with those of Ref. [45]. The transverse part of the p- meson 
self energy due to the nucleon loop is given as: 



(k) — nf DjT (k) + u.2 D T (k) + n^ D T (k) 

where, 



(25) 



\k\ 2 



|£f - k% - 4(E* N ) 2 k 2 - 4p 2 \k\ 2 \ ]nA _ 4K 2 k Q E* N M5 



( 9 2 p nn \ (&\ ( k pNN M* n \ r 
{ n 2 J \\k\J \ 2M N J Jo 



( 9jjNn\ 



\ 4tt 2 I \\k\z ) \AM 2 N J Jo E 



2p\k\ 3 
pdp 

TP* 
N 

00 p 2 dp 



2p\k\* 

n N (p*,T)+n N (p*,T))\nA, 



N 

/,-' + l^(^v) 2 - Ap 2 \k\ 2 \ h) ( I/.' 2 //, ■ 



2p\k\" 



n N (p*,T)+n N (p*,T)) 
In B 



2p|£f 



(26) 



The real and the imaginary parts of the self-energy can be obtained after performing the 



analytic continuation in both cases as k — > E + is where E = \J M 2 + |A;| 2 , with M p being 
the invariant mass of the p meson in the medium. The functions n p ™ (k) and IT^ (k) 
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acquire an imaginary part when A or B are negative. This happens when the variable of 
integration p, lies in the interval, 

l -\E^jl-{AM* N 2 )/M 2 -\k\\ <p< \\E^jl-{AM* N 2 )/M 2 + 



(27) 

The real and the imaginary parts of the logarithmic function A and B are given as: 
\nA = \n\A\-mQ(M 2 p -AM* N 2 ) and InB = In \B\ + m&(M 2 - AM* N 2 ). 

(28) 



B. piTTT interaction 

With the pirn interaction the p meson self-energy is calculated by using an effective 
Lagrangian density describing a system of charged pions and p mesons [22] 



J pTTTT 



\ | | 2 -\m\ | $ | 2 -\ PlH ,fr + \m 2 p^ 



(29) 

where $ is the complex charged pion field, p^ u = d^p u — d vPlx is the p field strength, and 
Dfi = — ig P mrPn is the covariant derivative. The one loop contribution to the polarization 
tensor is, 



n„„ (k) 



ig. 



piTTT 



d A p (2p + k)(2p + k) u 



(27r) 4 (p 4- k) 2 — ml [p 2 — ml 



O- 2 f d l 



d 4 p g^u 



(2tt)V -ml' 



(30) 



Separating IP" for the T = and temperature dependent contributions, we obtain the 
contribution for vacuum to be 



Klc (*) = - k 2 g^) 1 (^) 



1 - 



Ami 



k 2 



In 



5m 



1-^-1. 



A; 2 



(31) 
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where D is the renormalization constant. The renormalization constant D is fixed by 



Rell vac (k 2 = m 2 ) = in the free space. This gives D to be 



D = -2 



2 [ — ) In f^°±P° 



(32) 



where 2u = m p = 2ym 2 + p^. 

The T > contribution can be split into longitudinal and transverse parts as, 
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k 2 \ r°° p 2 dp 
47r 2 J \ \k 2 \J Jo UJ 
00 p 2 dp 



N(u) 



4- 



'4u 2 + fcg 
, 2p|fc| 



In a — 



/ 2/c a/ 
V p\k\ , 



In 6 



V 



-N(u) 



47T 2 / JO id 

2(k 2 + \k\ 2 ) (4p 2 + 2k 2 - \k\ 2 )\k\ 2 - k 2 (k 2 + 4oo 2 ) 



\k\ 2 



4p\k\-- 



ln a 



' k k 2 u s 
v p\k\ 3 , 



In 6 



(33) 



where, 



u = \Jp 2 + m 2 and iV(cu) 



Here is the pion mass. In the equation (31), a and b are defined by 

(k 2 + 2p\k\) 2 - Au 2 k 2 



(34) 



(k 2 -2p\k\) 2 -4u 2 ki 
k 4 - 4 (p\k\ + cuk 



2 ' 



k 4 -4{p\k\ -cuk ) 

The real and the imaginary parts of the self-energy can be obtained as before 
cleon loop, after performing the analytic continuation as ko — > E+ie where E 
The p meson spectral function is obtained as 

ImU L(T) (k) 



(35) 

for the nu- 
M 2 + \k\ 2 . 



A 



L(T) 



(k) 



M 2 - (m 2 + ReYi L[T) (k 



+ 



Imll 



L(T) 



(k) 



2 > 



(36) 



with Ul(t) being the total longitudinal (transverse) self-energy of p meson. 
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V. RESULTS AND DISCUSSIONS 



In this section, we present the results of our calculations of the spectral function of p- 
meson in the hot nuclear matter. In our calculations of the nuclear properties in the QHD-I 
model, we have used the values for the hadron masses in vacuum and the coupling constants 
as given in Ref. [43]. These values are g 2 = 109.626, g 2 = 190.431, m s = 0.52 GeV, 
M N = 0.938 GeV, m v = 0.783 GeV. In the chiral SU(3) model the parameters used are, 
m n = 0.1396 GeV, m K = 0.498 GeV, m w = 0.783 GeV, f n = 0.0933 GeV, f K = 0.122 GeV, 
Co = 0.10656 GeV, k = 2.37, h = 1.4, k 2 = -5.55, k 3 = -2.64, 5 = 2/33, xo = 0.4027 
GeV, g a N = 10.6, g^N = —0.47, g± = 2.7 [39]. The variation of the effective nucleon mass 
and effective chemical potential determine the spectral function of the rho meson in the hot 
and dense matter. The nuclear matter saturation density is chosen to be p — 0.16/m -3 [45] 
. The effective nucleon mass and chemical potential are plotted as functions of temperature 
for densities p# = 0,po,2p ,4p within the frameworks of QHD-I and chiral SU(3) in the 
mean field approximation. The results are shown in figures 1 and 2. Within QHD-I, for 
p B = as shown in figure la, is observed to remain almost a constant till T = 0.15 GeV 
above which it is seen to drop to around 0.82 at T = 0.2 GeV. But within chiral SU(3) 
model, as shown in figure 2a, is observed to drop to around 0.8 M N at a temperature, 
T = 0.18 GeV. A comparison of figures lb and 2b shows that the chiral SU(3) predicts 
higher values for the nucleon masses and chemical potential as compared to these in the 
QHD-I model. In figure 1, with the increase in density is seen to decrease for a given 
temperature. For a given density, is seen to increase with T till T = 0.15 GeV above 
which it decreases with temperature. For a given density, p* is observed to decrease with 
increase in temperature. Figure 2 follows the same trend. The results show that and p* 
in the RHA in both the models are seen to be quite similar to those obtained in the MFA, 
except that the values with RHA, are seen to be higher than the values of and p* with 
MFA. 

We consider isospin symmetric case, for which the longitudinal and transverse spectral 
functions are almost the same for all momenta. So for the sake of convenience, we show here 
only the longitudinal part of the spectral function. The parameters chosen in the calculation 
of the spectral function are, g 2 /An = 2.91, g 2 NN = 6.96, k p = 6.1. The tensor coupling is 
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FIG. 1: (Color online) Effective nucleon mass and chemical potential pf in MFA for QHD-I 
model in (a) for ps = 0, (b) for ps = po, (c) for ps = 2po, (d) for ps = 4po in GeV plotted as 
functions of temperature T in GeV. 

very important for the coupling of p meson with the nucleons. The coupling constants g P NN 
and k p are determined from the fitting to the nucleon-nucleon scattering data done by the 
Bonn group [43]. 

The longitudinal spectral function is studied for a particular p momentum ( \k\ = 0.75 
GeV ) for pb = and ps = po at T = 0.15 GeV in both the QHD-I and chiral SU(3) models. 
The results are shown in figure 3. In figure 3a, for ps = 0, the pure temperature effect (only 
pion loop) is seen to almost coincide with the mean field result in the Walecka model with 
the peak around 0.8GeV. With the nucleon loop alone, for small densities is large and 
the condition M 2 p > 4(M^) 2 is not satisfied. Hence there is no imaginary part, which means 
that there is only a real part to the spectral function implying that the spectral function is 
a 5 function peaked at M p = 0.77 GeV. When we consider both the loops within the RHA, 
it can be seen that the spectral function is seen to shift towards the low mass region and 
becomes a little sharper due to Dirac sea polarization. The p width decreases because of 
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FIG. 2: (Color online) Effective nucleon mass and chemical potential p* in MFA for Chiral 
SU(3) model in (a) for ps = 0, (b) for ps = po, (c) for ps = 2po, (d) for ps = ^po in GeV plotted 
as functions of temperature T in GeV. 

the decrease in the p mass. Figure 3c for the chiral model, is identical to figure 3a. Here 
though within MFA the peak is around 0.8GeV, within the RHA, the chiral SU(3) predicts 
slightly higher value for m* as a result of higher values. But for ps = po (figures 3b and 
3d), the RHA results for both the models give a spectral function which is further shifted 
towards the low mass region and is sharper. Here also the chiral model (figure 3d) predicts 
slightly higher values for m* and the p width. Within the MFA in both models the peaks 
are shifted to the high invariant mass region and becomes more wide since m* increases. 
This is consistent with the result of Ref. [38] for \k\ = 0. For smaller densities the RHA 
results are in agreement with the results in Ref. [45]. 

Consider the case when ps = 0. With the RHA in both the models (figures 4a and 4c) it 
can be seen that even though for T = 0.05GeV and T = O.lGeV the peaks coincide at around 
0.77GeV it shifts to around 0.7GeV and 0.72GeV for T = 0.15GeV as can be seen in figures 
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FIG. 3: (Color online) Spectral function against the invariant mass M p for \k\ = 0.75 GeV and 
T = 0.15 GeV for (a) ps = with the pion loop only, the pion and the nucleon loops within 
the RHA in Walecka model, the pion and the nucleon loops within the MFA in Walecka model, 
(b) pb = Po with the pion loop only, the pion and the nucleon loops within the RHA in Walecka 
model, the pion and the nucleon loops within the MFA in Walecka model, (c) ps = with the pion 
loop only, the pion and the nucleon loops within the RHA in chiral SU(3) model, the pion and the 
nucleon loops within the MFA in chiral SU(3) model, and (d) ps = Po with the pion loop only, the 
pion and the nucleon loops within the RHA in chiral SU(3) model, the pion and the nucleon loops 
within the MFA in chiral SU(3) model. The coupling constants are g p NN = 6.96 and k p = 6.1. 

4a and 4b. Hence there is larger reduction in the p mass when the is calculated with 
the Walecka model, as compared to in the chiral SU(3) model. Also with RHA the peaks 
become narrower since the p mass decreases. But with MFA as the temperature increases 
the peaks take slightly higher M p values in both the models implying a very slight increase 
in mass and a corresponding increase in the width. 

Keeping the 3- momentum of p meson fixed at 0.75GeV, the spectral function is plotted 
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FIG. 4: (Color online) Spectral function against the invariant mass M p for ps = 0, \k\ = 0.75 GeV 
for different values of T with Mjy and p* calculated with (a)RHA in Walecka, (b) MFA in Walecka, 
(c) RHA in chiral SU(3) and (d) MFA in chiral SU(3)models respectively. The coupling constants 
are g p NN = 6.96 and k p = 6.1. 

for various T values for a higher value of the baryon density ps = Po m figure 5. Here again 
figures 5a and 5c are almost identical, figure 5c which corresponds to the RHA in chiral 
SU(3) model gives slightly higher values for m* and the p width. In both the models within 
RHA, m* and p width increase slightly with temperature. This is consistent with the result 
in Ref . [38] . But in figure 5b with MFA and Mjy calculated in the Walecka model the spectral 
function is observed to have double peaks at T = and at T = 0.05 GeV. But at a higher 
temperature the spectral function becomes a smooth curve. The p width is a maximum 
at T = 0.1 GeV. Using MFA, in the chiral SU(3) model, the peak of the spectral function 
(shown in figure 5d) is almost at around the same M p as in Walecka model shown in figure 
5b, although in the high invariant mass region. The p width is observed to be maximum at 
T = 0.05 GeV in figure 5d. The spectral function becomes much wider in comparison with 
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FIG. 5: (Color online) Spectral function against the invariant mass M p for p B = po, \k\ = 0.75 
GeV for different values of T with Afjy and p* calculated with (a)RHA in Walecka, (b) MFA in 
Walecka, (c) RHA in chiral SU(3) and (d) MFA in chiral SU(3)models respectively. The coupling 
constants are g p NN = 6.96 and n p = 6.1. 

the RHA result. 

For Pb = 2p , the result for the spectral function is shown in figure 6. In figures 6a 
and 6c with the RHA calculation in the Walecka and chiral models respectively, the plots 
are observed to be very similar. But a comparison with figures 5a and 5c, which are for a 
lower density, p B = Po, shows that in figures 6a and 6c the peak is at a lower M p value for 
all temperatures. In figure 6a, with the increase in temperature the peaks almost coincide 
with the peak at T = 0. In figure 6c there is a very slight increase in m* with increase 
in temperature. With MFA in the Walecka model (figure 6b) it can be seen that as the 
density is increased to p B = 2p the double peaks as seen for p B = p in figure 5b, merge to 
a single peak, which is further shifted to the high mass region and the width increases with 
increasing T. This suggests an increase in m* and a corresponding increase in the width. 
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FIG. 6: (Color online)Spectral function against the invariant mass M p for ps = 2po, \k\ = 0.75 
GeV for different values of T with Afjy and p* calculated with (a)RHA in Walecka, (b) MFA in 
Walecka, (c) RHA in chiral SU(3), and (d) MFA in chiral SU(3) models respectively. The coupling 
constants are g p NN = 6.96 and n p = 6.1. 

In figure 6d with MFA in the chiral model, there is a further shift in m* to the high mass 
region due to a higher value of M^. 

For a higher density of ps = 4po, the result for the spectral function is shown in figure 
7. A comparison with figure 6 shows that within the RHA in both the models the spectral 
function peaks at T = shift slightly to higher M p values but, with the increase in tem- 
perature the peaks almost coincide at a slightly lower value approaching a 5- like function. 
A comparison with figures 5 and 6 shows that within RHA, a reduction of the p mass is 
found upto around 2p . At higher densities, the density dependent part of the p meson 
self energy, describing the Fermi sea fluctuations starts to be more dominating, leading to 
slightly increasing masses. Within MFA in the Walecka model (figure 7b), it can be seen that 
the contribution is small. Since is very small, at ps = 4po, there will be an imaginary 
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FIG. 7: (Color online)Spectral function against the invariant mass M p for ps = ^po, \k\ = 0.75 
GeV for different values of T with Afjy and p* calculated with (a)RHA in Walecka, (b) MFA in 
Walecka, (c) RHA in chiral SU(3) and (d) MFA in chiral SU(3)models respectively. The coupling 
constants are g p NN = 6.96 and n p = 6.1. 

part to the p self energy due to the nucleon loop. The inclusion of this reduces the spec- 
tral function in our region of interest. But in figure 7d, within MFA in the chiral SU(3) 
model, Mjy is still too large to give an imaginary part. So again the pion loop alone gives 
an imaginary part which is responsible for the spectral function in figure 7d. 

The longitudinal spectral function is studied keeping the density at ps = for different 
temperatures by changing the p momentum \k\ to 0.25GeV in figure 8. With the RHA in 
both the models (figures 8a and 8c) it can be seen that even though for T = 0.05GeV and 
T = O.lGeV the peaks coincide at around 0.77GeV, it shifts to a value around 0.72GeV 
in figure 8a and around 0.75GeV in figure 8c for T = 0.15GeV. Comparison with figures 
4a and 4c suggests that when the 3-momentum of the p is decreased the mass drop is 
smaller. But with MFA, for \k\= 0.25GeV, as the temperature increases, the peaks take 
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FIG. 8: (Color online) Spectral function against the invariant mass M p for ps = 0, \k\ = 0.25 GeV 
for different values of T with M*^ and p* calculated with (a)RHA in Walecka, (b) MFA in Walecka, 
(c) RHA in chiral SU(3) and (d) MFA in chiral SU(3)models respectively. The coupling constants 
are g p NN = 6.96 and k p = 6.1. 

higher M p values in both the models implying a slight increase in mass and a corresponding 
increase in the width. Within MFA in both models the finite temperature effects increase 
with decreasing momentum while with RHA the finite temperature effects decrease with 
decreasing momentum. 

For density, ps = po, the p momentum fixed at \k\ = 0.25GeV, the spectral function is 
plotted for different temperatures in figure 9. Figure 9a which corresponds to the Walecka 
model within the RHA gives slightly higher values for the peak position in comparison to 
the values obtained with \ k\ = 0.75GeV as shown in figure 5a. The plots for T = 0, T = 0.05 
GeV, and T = 0.1 GeV almost coincide, but the plot for T = 0.15 GeV shifts to higher 
M p value and the peak becomes broader suggesting an increase in m* and the p width. A 
comparison with figure 8a shows that for the same p 3-momentum and temperature, m* and 
p width decrease with increasing value of ps- This significant reduction of the p mass due 
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FIG. 9: (Color online) Spectral function against the invariant mass M p for ps = po, \k\ = 0.25 
GeV for different values of T with Afjy and p* calculated with (a)RHA in Walecka, (b) MFA in 
Walecka, (c) RHA in chiral SU(3) and (d) MFA in chiral SU(3)models respectively. The coupling 
constants are g p NN = 6.96 and k p = 6.1. 

to Dirac sea polarization is found upto around nuclear saturation density. Moreover for the 
same p 3- momentum, if the density is increased from Pb = to Pb = Po, m * p and p width 
slightly increase with increase in temperature in the latter case, but are observed to decrease 
with temperature in the former case. Figure 9c corresponding to the chiral SU(3) within the 
RHA, gives higher values for both m* and the width than as compared to Walecka model 
shown in figure 9a. The maximum value of m* and width is for T = 0.15 GeV, the value of 
m* being around 0.52GeV. Figure 9b and 9d correspond to the Walecka and chiral SU(3) 
models within the MFA. In MFA, m* and p width tend to increase very much since we are 
considering only the Fermi sea fluctuations. Our results are similar with the result in Ref. 
[38]. Also a comparison with figures 5b and 5d suggests that finite temperature effects are 
remarkable at low values of 3- momentum of p -mesons in MFA. 

With the 3- momentum of p- meson as 0.25GeV, for ps = 2po, the spectral function 
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is plotted for different values of temperature T in figure 10. Comparing figures 10a and 
9a corresponding to Walecka model within RHA, it can be seen that at higher density 
the spectral function follows a different pattern. For T = the peak is around 0.62GeV. 
With the increase in temperature the peak shifts further towards low mass region and the 
width also decreases correspondingly contrary to what is observed in figure 9a. In RHA, 
a significant reduction of the p mass due to Dirac sea polarization is found up to around 
nuclear saturation density. But at higher densities, the density dependent part YI D41U of the 
vector meson self energy becomes more dominant, leading to an increase in m*. Comparison 
of figures 10a and 6a shows that for the same ps, m* and the p width increase with decrease 
in 3-momentum of the p meson. Figure 10c for the chiral SU(3) model within RHA also 
follows the same pattern as figure 10a with the peaks taking slightly higher M p values. 
Figures 10b and lOd correspond to the Walecka and the chiral SU(3) models within the 
MFA. The spectral function contribution is very small in figures 10b and lOd in the region 
of temperatures we have considered. The peak may be further shifted towards still higher 
invariant mass region indicating a very high value of m*, outside the range of M p plotted 
here. Comparison of figures 10b and 6b indicates an increase in the p mass in figure 10b as 
compared to figure 6b in the Walecka model. The same happens for figures lOd and 6d also 
for the case of chiral SU(3) model. 

With the nucleon density, pb = 4p , keeping the p momentum fixed at 0.25GeV, the 
spectral function is plotted for different temperatures in figure 11. In figures 11a and 11c 
corresponding to the Walecka and the chiral SU(3) models within RHA, the results are very 
similar. The spectral function peaks shift to higher M p values and the width increases with 
increasing temperature. Within MFA in the Walecka and the chiral SU(3) models (figure 
lib and lid) it can be seen that there is no contribution to the spectral function in the 
region of our interest. Probably the peak may be further shifted towards the high invariant 
mass region than in figures 10b and lOd. 

So far we have studied the spectral function of the p meson using the p — N coupling 
strengths as determined from the NN forward scattering data. Now we study the p meson 
spectral function in the chiral SU(3) model with the nucleon-p coupling, g P NN, as determined 
from the symmetry relations [46]. Here we take the tensor coupling as a parameter in our 
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FIG. 10: (Color online) Spectral function against the invariant mass M p for ps = 2po, \k\ = 0.25 
GeV for different values of T with Mjy and p* calculated with (a)RHA in Walecka, (b) MFA in 
Walecka, (c) RHA in chiral SU(3) and (d) MFA in chiral SU(3)models respectively. The coupling 
constants are g p NN = 6.96 and n p = 6.1. 

calculations since this coupling cannot be fixed from infinite nuclear properties. The coupling 
constants chosen are q p nn = 4.27 and k p nn = 2,4,6 respectively. Here first we consider 
the case when \k\ = 0.75GeV and Pb = Po and study the spectral function by varying the 
temperature. The results are shown in figure 12. Figures 12a, 12c and 12e correspond to the 
the chiral SU(3) model within RHA with the tensor coupling constants as 2,4,6 respectively. 
In figure 12a where the vector coupling is larger than the tensor coupling we can see that the 
peak is around 0.5GeV. When the temperature increases the peak slightly shifts indicating 
a very small increase in m* and p width. In figure 12c, the tensor and the vector couplings 
are comparable. Here it can be seen that the peak is further shifted to the low invariant 
mass region to a value around 0.4GeV and the behaviour is seen to be similar to figure 12a. 
Figure 12e correspond to the case where the tensor coupling is dominating over the vector 
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FIG. 11: (Color online) Spectral function against the invariant mass M p for ps = 4poj |^| = 0.25 
GeV for different values of T with Mjy and p* calculated with (a)RHA in Walecka, (b) MFA in 
Walecka, (c) RHA in chiral SU(3) and (d) MFA in chiral SU(3)models respectively. The coupling 
constants are g p NN = 6.96 and n p = 6.1. 

coupling. The peak shifts much towards the low mass region around 0.33GeV and becomes 
very narrow. Figures 12b, 12dand 12f correspond to the MFA case. In figure 12b, where the 
vector coupling is dominant almost all the peaks coincide with the peak around 0.82 GeV 
and the p width is also slightly reduced. In figure 12d, the vector and the tensor coupling 
are comparable. Here we can see that as the temperature is increased the peak is shifted to 
the high invariant mass region. The p width is maximum at T = O.lGeV above which it is 
further reduced. In figure 12f, where the tensor coupling is dominant the spectral function 
gets shifted with temperature again to the high mass region and the width increases with 
temperature. It can be concluded that within MFA, the p is more stable when the vector 
coupling is more dominant. But within RHA, the p meson is more stable when the tensor 
coupling is predominant. 
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FIG. 12: (Color online) Spectral function against the invariant mass M p for ps = po, \k\ = 0.75 
GeV for different values of T with (a)RHA with h p n n = 2, (b) MFA with k p nn = 2, (c) RHA with 
KpNN = 4, (d) MFA with k p nn = 4, (e)RHA with k p nn = 6, (f) MFA with k p nn = 6 respectively. 
9 P NN = 4.27. 

We then study the spectral function varying ps keeping the \k\ = 0.75 GeV and T = 
0.15GeV. The results are given in figure 13. Figures 13a, 13c, and 13e correspond to the 
RHA calculation. With figure 13a (k p nn = 2 ) the peak shifts towards the low mass region 
as pb increases and the width is also reduced. In figure 13c (k p atv = 4) the reduction in 
m* with is more and the peak is a 6 function. The p width is also reduced as compared 
to figure 13a. Figure 13e (k p nn — 6) shows that m* and p width are reduced very much 
with increasing p#. Figures 13b, 13d and 13f correspond to calculations in MFA. In figure 
13b, as p# increases the peak is seen to be shifting to the high invariant mass region. In 
figure 13d (k p nn = 4) the shift and the width are observed to be larger. For higher densities 
the contributions are seen to be negligible. With k p nn = 6 the peak is still in the high 
invariant mass region for Pb = Po and the contributions are negligible for higher densities in 
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FIG. 13: (Color online) Spectral function against the invariant mass M p for T = 0.15GeV , \k\ = 
0.75 GeV for different values of the nucleon density ps (a)RHA with k p nn = 2, (b) MFA with 
KpNN = 2, (c) RHA with k p nn = 4, (d) MFA with K pNN = 4, (e)RHA with K pNN = 6, (f) MFA 
with k p nn = 6 respectively. g p NN = 4.27. 

our region of interest. Within RHA, m* and the p width decrease with increasing nuclear 
density whereas within MFA both increase with increasing nuclear density. 

The spectral function is also studied by varying \k\, keeping the temperature and ps fixed. 
The results are given in figure 14 for T = Q.lbGeV and ps = Po- When k p nn = 2 with 
RHA (figure 14a) the spectral function has a peak at around 0.62GeV for \k\ = 0.25GeV. As 
the momentum increases the peak shifts and almost remains at the same position at around 
0.52GeV. The p width also decreases. With increasing tensor couplings (figure 14c and 14e) 
the peaks shift further towards the low invariant mass region and become narrower. With 
MFA, when k p nn — 2, the peak shifts to higher values upto \k\ = 0.5GeV after which it 
again falls to around 0.8GeV. When h p nn = 4 m* increases upto \k\ = 0.75GeV and at \k\ 
= lGeV the peak is around 0.85GeV and the width is reduced. With k p nn = 6 the shift 
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FIG. 14: (Color online) Spectral function against the invariant mass M p for T = 0.15GeV, and 
Pb = Po for different values of \k\ (a)RHA with k p nn = 2, (b) MFA with k p nn = 2, (c) RHA with 
KpNN = 4, (d) MFA with k p nn = 4, (e)RHA with k p nn = 6, (f) MFA with k p nn = 6 respectively. 
9 P NN = 4.27. 

and the width increase with \k\ indicating that the particle is highly unstable. 
VI. SUMMARY 

To summarize, in the present work, we have investigated the temperature and density 
effects on the p meson spectral function with the effective lagrangian in the ambit of the 
QHD-I model and the chiral SU(3) model in isospin symmetric pion and nucleon media. 
With only the pure temperature effect, i.e., with pion loop only, the medium corrections 
are observed to be modest even upto a temperature, T = 0.15 GeV. But the inclusion 
of the nucleon loop drastically changes the medium properties of the p meson, meaning 
that the density effect is a remarkable effect as compared to the pure temperature effect. 
In QHD-I model, has a larger drop with density as compared to the chiral SU(3) 
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model. This is reflected in the spectral function of the p meson. Even for small densities 
the spectral function is shifted towards the low invariant mass region significantly and the 
spectral function becomes very sharp in the RHA. For the same density the spectral function 
gets shifted slightly towards the high invariant mass region and becomes broad in the MFA. 
This is consistent with the recent experimental result which supports a broadening induced 
in the medium. 

Also it can be seen that the spectral function depends on the coupling strengths of the p 
meson with the nucleons. If the tensor coupling is dominant then m* has a tendency to take 
comparatively lower values. If the vector coupling is predominant then m* has a tendency 
to take slightly higher values in RHA. The p width also follows the same pattern. But in 
the MFA when the vector coupling is dominant m* and p width take comparatively lower 
values than when the tensor coupling is dominant. 

— * 

The spectral function of the p meson also depends on its momentum \k\. For the same 
temperature and baryon density, ps, the spectral function shifts to the low mass region 
indicating a drop in m* with \k\ in RHA irrespective of whether the vector coupling or the 
tensor coupling is dominant. In MFA for lower momenta the peak is at a higher M p value, 
but as the momentum increases the peak takes a comparatively lower value though in the 
high invariant mass region . 
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